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The real-space RG approach is applied to study critical temperatures of system 
consisting of interacting spin chains of spin 5=1 with an inner antiferromagnet- 
ic exchange which form a honeycomb crystal lattice. Using anisotropic Heisenberg 
model we calculate critical temperature as a function of anisotropic parameter and 
the ratio of interchain and intrachain interactions. A comparison our results with 
those obtained from RGRS calculations for the same model of spin-1/2 on a square 
lattice is given. 



I. INTRODUCTION. 



The study of low-dimensional magnetic systems made of spin chains attracts much at- 
tention for last decades. If the interchain coupling is less than intrachain one, at higher 
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temperatures these systems exhibit properties intrinsic to one-dimensional magnets. How- 
ever, at lower temperatures the interchain interaction comes into play and begins to govern 
a magnetic behavior of the system. 

In this work we apply the real space renormalization group (RSRG) procedure to the 
antiferromagnetic Heisenberg (AH) model of spin S = 1 on a honeycomb lattice made of 
antiferromagnetic spin S = 1 chains with an antiferromagnetic interchain coupling. We 
are interested in the critical properties of the 2D system, and explore how the interchain 
coupling affects the criticality of the model. Simultaneously, we elucidate an effect of an ex- 
change anisotropy on the critical temperature and build the phase diagram of the anisotropic 
Heisenberg (AH) model of spin S = 1 on a honeycomb lattice with two kind of exchange 
couplings. 

First of all, the interest to the problem is motivated by synthesis of a family of related 
organic biradicals PNNNO, PIMNO and F 2 PNNNO whose magnetic properties have 
been examined by susceptibility and magnetization measurements^. Each biradical involves 
two spins of S — 1/2, which are coupled ferromagnetically (Jf)- In their turn, these spin pairs 
couple antiferromagnetically in the crystal. Due to the strong ferromagnetic coupling Jf, 
antiferromagnetic chains of S = 1 are seen in both PNNNO and F2PNNNO (see Fig. 2(a) 
Type-I in Ref.— ). PNNNO has interchain interactions in three-dimensions (3-D), whereas 
F 2 PNNNO has two-dimensional (2D) interchain interactions. PNNNO is well understood 
by one-dimensional antiferromagnetic chain model. The compound undergoes 3D Neel order 
at ~ 1 K due to the weak interchain coupling. The crystal structures of F 2 PNNNO is 
shown in FigQJ In the extreme limit of Jp — > 00, the model becomes equivalent to the 
coupled antiferromagnetically (J'af) antiferromagnetic uniform chains of S — 1 with the 
intrachain exchange integral Jaf- F 2 PNNNO with comparable values of two kinds of the 
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antiferromagnetic interactions Jaf and J' AF forms a 2D system on a honeycomb lattice. It 
has been found that a magnetism of F 2 PNNNO is characterized by the singlet ground state 
and an energy gap above the state. This finding is supported by the high-field magnetization 
process which shows a plateau in the magnetization curve. We have to point out once that 
a modelling of real spin- 1/2 system by spin 1 system implies some mistake and one have to 
apply the results of the RSRG analysis to the real compounds with caution. 

The results for the spin-1 lattice may be useful also for the theory of S = 1 bosons 
(ultracold atoms 23 Na) trapped in an optical lattice in the regime of one particle per site 
for suitable interaction between the bosons^. 

According to recent quantum Monte Carlo results for spin- 1/2 weakly anisotropic an- 
tiferromagnets on the square lattice an ordered low-temperature phase develops for very 
small anisotropy of order 1CT 3 -r- 10~ 2 (in units of exchange coupling)^. These results are in 
contradiction with the RSRG treatment for the same model predicting considerably larger 
value of the critical anisotropy (~ 0.2). In view of this discrepancy, the RSRG analysis for 
the S > 1/2 case is of a particular theoretical interest and besides it is instructive to com- 
pare a critical behavior of systems of integer and half-integer spins. We mention the results 
of investigations of ground states of spin-3/2 (spin-2) systems on the hexagonal (square) 
lattice by using vertex state models^. It has been shown that these systems exhibit an 
anisotropy-induced second order phase transition from disordered phase with exponential 
decay of two-spin correlations to the Neel order phase with dominant long-range correlations. 
The phase transition belongs to the 2D Ising-model universality class. 

After the RSRG approach has been applied successfully to study the 2D Ising systems^ 
a number of works have been dedicated to the investigation of phase transitions in quantum 
systems within the method^^i. In last decade, the RSRG methods have been performed to 
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calculate the phase diagram for the anisotropic antiferromagnetic Heisenberg (AAH) model 
of S = 1/2 on the square lattic o 12 ! 13 . This approach uses a hierarchical lattice to approximate 
the square one, and performs a partial trace over internal degrees of freedom. Recently, in 
order to study weakly interacting classical and quantum spin chains the linear-perturbation 
real-space renormalization group (LPRG) method has been suggested^. The LPRG uses 
the existence of the small parameter - the ratio of interchain to intrachain coupling. This 
perturbative method involves RG transformation, which for the Ising spins is the standard 
decimation procedure and for the quantum spins is the generalization of the Suzuki-Takano 
approximate decimation^. However, in practice, the LPRG using perturbation theory with 
the interchain interactions as the perturbation parameters is only reliable for small values 
of the ratio interchain to intrachain interaction. 

We use an extension to S = 1 case of the quantum RSRG approach originally suggested 
by Mariz et al. for the 5=1/2 AH model^. An application of the renormalization group 
method to interacting quantum spin chains encounters standard difficulties connected with 
a necessity of a decomposition of the exponential operators, and additional proliferation of 
the new interactions due to the vector character of the spin operators. 

The paper is organized as follows. In Sec. II the RSRG method is formulated for S = 1 
and applied to the honeycomb lattice formed by interacting antiferromagnetic chains. Our 
results and conclusions are given in Sec. III. 
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II. MODEL. 

We consider a system (anisotropic Heisenberg model) whose dimensionless Hamiltonian 
is defined by 

-(3H= £ K t3 

<i,j> 

where (3 = l/fc^T, K%j = Jij/ksT, Jij is the exchange coupling constant, (ij) denotes first- 
neighboring lattice sites, Ay is the anisotropic parameter, and S" {a = x, y, z] is the spin 
1 on the site i. The Hamiltonian jH) describes the Ising (Ay = 1), isotropic Heisenberg 
(Ay = 0) and XY-model (Ay = -oo). 

We mention briefly main operations lying in base of the approach of Ref.— . A parallel 
array of two bonds characterized respectively by (Ki,Ax) and (K 2 ,A 2 ) is equivalent to a 
single bond characterized by (K p , A p ) given by 

K P = K 1 + K 2 , 

K P A P = KxAx + K 2 A 2 . 

The extension of the approach, that is nothing but Migdal-Kadanoff procedure^ 1 ^, to n 
parallel bonds is straightforward. 

For a series array of two bonds characterized by (Ki, A x ) and (K 2 ,A 2 ) the procedure 
involves difficulties due to non-commutativity effects. Under rescaling and removal of inter- 
mediate spins (decimation), the Hamiltonian changes and is characterized by a new set of 
parameters that are functions of the original set. The initial Hamiltonian is given by 

H 123 = K x [(1 - A x ) + SfS!) + SfSI) + K 2 [(1 - A 2 ) (S£Sf + Sf Sf) + SfSf] . 

We have to replace this array by a single bond whose Hamiltonian is 

H' 12 = K s [(1 - AO (5?5J + S!S y 2 ) + SIS*} + K' . (2) 



1 - Ay) (SfSf + SfS]) + SIS* 



(1) 
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We impose the preservation of the partition function between terminal sites 1 and 2, i.e. 

exp H[ 2 = Tr 3 exp H 132 , (3) 

where Tr 3 denotes the tracing operation over the states of the intermediate spin 3, K' is an 
additive constant included to make Eq. Q possible. Equation Q establishes the relation 
between the set of parameters (K\, Ax), (K 2 , A 2 ) and the set of renormalized parameters 
(K s , A s , K' ). For the anisotropic spin-1/2 Heisenberg model corresponding expressions may 
be found explicitly^,—. 

In order to construct recursion relations of the RG transformation Q for the spin-1 model 
we should expand both sides over an appropriate matrix basis to equate coefficients of the 
expansion. To avoid a proliferation problem we chose those quantities that correspond to 
coupling constants in the initial Hamiltonian JH). 

Due to properties of Pauli matrices the matrix representation for H and exp H in the case 
of spin 1/2 has the same form, i.e. H and exp if have non-zero matrix elements in the same 
positions. For the 5 = 1 this rule does not hold that results in an essential complication of 
RG procedure. 

Following the treatment of Ref we expand exp H[ 2 as 

oo oo 

ex P #( 2 = £ J2 K nin2 A^®A n 2 \ 

711=0 T12=0 

where £g> is the outer product, v4 12 is ordinary powers of the spin operators S^f' 2 and 
the coefficients K nin , 2 depend on K S ,A S and K' . Since, A\ 2 are the 3x3 matrices we 
expand them over the basis that consists of the polarization matrices T g fc (k = 0,1,2 and 
q = —k, —k + 1, . . . , k) (see Appendix A) 

A? = a{T$+ £ b M (T i ) 1 M + £ c M (T&, (< = 1,2). 

M=±1,0 Af=±2,±l,0 



In their turn, the matrices T may be written through the spin operators explicitly^* 



T ° = 4^ T li = =4 ( S * ± iS v) > T c 



v/3 



T 2 ±l = T \ \{S x S z + S Z S X ) ± i (S y S z + S z S y )} ,T 2 ±2 = \ [(S x ) 2 - (S y ) 2 ± z (S^ + S y S x ) 

The transformation exp (if) should preserve the symmetry of the Hamiltonian H. Thus, the 
requirement of invariance gives the most general form 

exptf 12 = «i ((TO? ® (T 2 )J)+a 2 ((1$ ® (T 2 )J)+a 3 ((TO J ® (T 2 )j!)+/3 ((T)J ® (T 2 )^ + (T\) 
+ 7 ((Ti) 2 ® (T 2 ) 2 2 + (T x ) 2 _ 2 ® (T 2 ) 2 ) + a ((T) 2 ® (T 2 ) 2 x + (T) 2 ^ ® (T 2 ) 2 ) (4) 

with a new set of coupling parameters «i, a 2 , a 3 , /3, 7, a. For the Hamiltonian (j2j) the above 
formula gives 



exp H 



12 



\ 



A 12 
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are 


A 12 = 


E «i/3 + 


a 2 /2 + a 3 /6, 5 12 


= a x /3 - a 3 /3, C 12 = 



-13/2 - a/2, D 12 = at/3 - a 2 /2 + a 3 /6, E 12 = at/3 + 2a 3 /3, F 12 = -(3/2 + a/2, G 12 = 7. 

Similarly, we obtain the closed expression for the expansion of Tr 3 exp ifi 23 akin to that 
of exp H12 with the coefficients in the expansion (j3J are functions of the parameters coming 
into f/i 23 (see Appendix B for details). 
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To calculate the exponentials we will diagonalize numerically the 9x9 and 27 x 27 
matrices associated with Hi 2 and H123. By using 

expH 12 = ^12 exp(i^) 6^2, exp#i 23 = U 123 exp(H? 23 )Ul 23 , 

where Ui 2 , U\ 23 are the unitary matrices turning Hi 2 , Hi 23 into the diagonal forms H® 2) H® 23 , 
we can find numerically exp H\ 2 and exp ifi 23 as functions of corresponding coupling param- 
eters. The same matrix structure © of exp H\ 2 and Tr 3 exp Hi 23 is supported by the nu- 
merical calculation. The numerical procedure produces a set {oti, a 2 , ct 3 , (3, 7, a} for exp H\ 2 
and \ol\i a 2 , a 3 , (3, 7, a\ for Tr 3 exp if 132 . To obtain the required RG equations we impose 

ai = A 12 + B l2 + D12 = A12 + Bn + D12 = a h (6) 

a 2 = A 12 - D 12 = A l2 - D 12 = a 2 , (7) 

(3 = -C 12 - F 12 = -C 12 -F 12 = p (8) 

and 

«3 = E\ 2 — B\ 2 = E\ 2 — B12 = «3, 7 = G12 = G12 = 7, o" = F12 — C*i 2 = F12 — C12 = o". (9) 

The number of these equations exceeds the number of interactions that enter into the Hamil- 
tonian J2J) because all possible bilinear couplings between terminal sites come into play. 
Thus, in order to carry out the RG decimation we retain the three equations (|6I7I8|) which 
implicitly determine K s , A s and K' as functions of (Ki, A x ), (K 2 , A 2 ). This set of equations 
is a counterpart of RG relations for the case of S = 1/2 (see Eqs.(12) in Ref.-^). 

Now, we have to chose an appropriate hierarchical lattice. We take one of the simplest 
cells, conserving a point symmetry of the full lattice, with 6 sites and 6 bonds, depicted in 
FigI2J We then impose that the correlation function between the two terminal sites 3 and 6 
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of the original and renormalized graphs are preserved. At the first step we apply decimation, 
the spins 1 and 3 (or 4 and 6) survive whereas the spins 2 and 5 are removed. At the second 
step the decimation procedure is repeated removing the spins 1 and 4. Finally, to obtain the 
renormalized parameters we apply Migdal-Kadanoff bond moving combining the "pieces" in 
parallel, which leads to the recursion relations 

(K s ,A s ) = Rs (K 2 ,A 2 ;K 1 ,A l ), 

(k' s ,A' s )=R s (K s ,A s ;K 1 ,A 1 ), (10) 
(K P ,A P ) = 2(K' S ,A' S ). 
We have evaluated numerically the renormalization transformation from the original 
set of coupling parameters to the set of renormalized parameters. Critical points are then 
evaluated as non-trivial fixed points of the above relations which can be rewritten as the 
composite function 

(K p , A p ) = 2R S (R s (K u A i; K 2 , A 2 ) ; K u A,) . (11) 

Unlike to the case of S = 1/2 we can not obtain RG relations explicitly. Instead of 
this, we briefly outline the numerical procedure. Input fixed parameters are the ratio of 
interchain to intrachain coupling C\ = J'af/Jaf arid the anisotropic parameter A = A x 
and what's more A 2 = C\A (the feature of anisotropy is the same both for the intrachain 
and interchain couplings). At given starting value Ku of the intrachain coupling (then the 
interchain coupling is K 2 i = CiKu) we apply two successive decimation steps to produce a 
renormalized coupling K' s . During each of the step we solve Eqs. (lfil7l8j) using the standard 
routine for non-linear systems of equations^. Then, we double the result obtained after these 
transformations to get a final value Kj depending on the starting value Ku. To complete 
we find a fixed point K c of the equation Ku = Kf(Ku) by using bisection method. 
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III. RESULTS. 

The critical inverse temperature K c = 1/T C as a function of C\ for several A values 
is presented in FigO As seen, the critical temperature rapidly decreases when the inter- 
chain coupling becomes weaker. Our results for the critical temperature as a function of 
the anisotropy A are shown in FigHJ The universality class for the whole critical curve is 
the same as for the Ising model. By contrast to some RSRG calculations for S = 1/2 the 
phase diagrams for ferromagnetic and antiferromagnetic models are the same: the critical 
temperature reaches zero at a critical value of A, A c , which is greater than zero. The weaker 
interchain coupling the stronger quantum fluctuations. So, as one could expect, the A c value 
is larger if the C\ shifts to lower values. 

For the lowest temperatures we could work we have observed no sign of the reentrant 
behavior found in some previous RGRS treatments. The Neel temperature behaves as 

T 1 

± AT rs - / 

N ln(A-A c ) 

near A = A c that agrees with the result found for the case S = 1/2. Our calculation can not 
be carried out down to T = 0; therefore we can not make any definite conclusions about the 
ground state of the model. The scaling law holds for different values of the ratio of interchain 
to intrachain coupling (FigEJ). In Ref . 21 ' 22 , the logarithmic dependence of T/v and T c with 
respect to A — A c is established using scaling arguments with A c = 0. Recent quantum 
Monte Carlo results^ for the anisotropic 2D S = 1/2 Heisenberg model have shown that it 
develops an ordered low-temperature phase even for very small anisotropies A ~ 10~ 3 , 10 -2 . 
The latter gives strong evidence of large values of a critical anisotropy is an artifact of the 
real-space renormalization approach. 

At this point, it is worthwhile to compare our results with those from RGRS calculations 
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for antiferromagnetic AH model of 5 = 1/2 on a square lattice. These RSRG analysises lead 
to non-equivalence between the criticality of the ferromagnetic (F) and antiferromagnetic 
(AF) models, a reentrant behavior in the (T, A) diagra m 12 ' 1 ^ . 

It is well known that in classical spin models, such as as the Ising or classical Heisenberg 
models, on bipartite lattices the critical temperature is the same for ferromagnetic exchange 
(Curie temperature) as for antiferromagnetic exchange (Neel temperature). This is a direct 
consequence of the free energy being an even function of the exchange parameter. However, 
for the quantum spin 1/2 Heisenberg model the Curie and Neel temperatures are unequal. 
Recently, this question has been reinvestigated using high-temperature series expansions 
for the spin 1/2, 1 and 3/2 Heisenberg ferromagnet and antiferromagnet in 3-dimensions 2 ^. 
The difference between the temperatures decreases rapidly with increasing spin 5. In some 
quantum systems, such as the quantum spin 1/2 XY and transverse Ising models, an iso- 
morphism between the criticality of the F and AF cases is observed^. Critical properties 
of the quantum spin-1/2 2D Heisenberg model with anisotropic interaction treated by the 
Green's function technique yields T c = T N for all values of the anisotropy paramete r 26 ! 27 . 
RSRG methods give contradictory results of the problem because of underlying approxima- 
tions whose effects are hard to control in a systematic way. In Refs . 12 ' 1 ^ it was obtained 
by RG approach T N < T c for 2D anisotropic Heisenberg limit < A < 1 due to a special 
choice of the hierarchical lattice approximating a square one. The critical temperature T c 
for 2D ferromagnetic spin 1/2 anisotropic Heisenberg model tends gradually to zero when 
decreasing the anisotropy parameter A, i.e. T c = in the isotropic Heisenberg limit A = in 
accordance with Mermin- Wagner theorem 2 ^ The results for the antiferromagnetic exchange 
are very similar to the 5 = 1/2 case, there is no long-range Neel order for the anisotropy 
parameter A < A c . The question is not yet settled and more work is needed to put these 
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points on firmer grounds. 

An observation of reentrance behavior for the spin 1/2 antiferromagnetic AH model on 
the square lattice has been reported by some author o 12 ' 1 '? . This result suggests that there 
is an ordered phase at relatively high temperature but not at very low temperature. A full 
understanding of this phenomena is still lacking, but it's most likely the reentrance behavior 
is an artifact of the RGRS method. In Refs.— the authors attribute the reentrance to the 
effect of finite size in the renormalization and for larger clusters it should be absent. 

In conclusion, the real-space renormalization group is employed to study the anisotropic 
Heisenberg model of spin S = 1 on a honeycomb lattice with two kinds of antiferromagnetic 
couplings. We calculate dependencies of the critical temperature on the parameters of the 
magnetic anisotropy, and on the ratio of interchain and intrachain exchange interactions. 
The entire critical line is found to belong to the universality class of the Ising model. In 
accordance with the early RSRG predictions for the antiferromagnetic AH model of spin S = 
1/2 on a square lattice our calculations recover an existence of large finite critical anisotropy 
A c that should be considered, however, as an artifact of the real-space renormalization 
technique. 
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Appendix A 

The polarization matrices for the spin S = 1 has the explicit form 
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We first take care of H 12 and express it in the basis \MiM 2 ). In this basis H 12 becomes 
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where Ws — Kg (1 — As). 

Performing the same calculation for -^123, now using the basis |MiM 3 M 2 ), we arrive at 

a 27 x 27 which has 4 independent blocks of size 9x9 
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where the 9 x 9 A±, Bi, B 2 , A 2 and A 3 are given by 
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B 2 = B l 
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Wi 

where W l = K x (1 - A x ) and W 2 = K 2 (1 - A 2 ). 
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Phc. 1: The magnetic model for F2PNNNO. (a) Uniform chains with intramolecular ferromagnetic 
coupling ( Jf) and intrachain antiferromagnetic coupling (Jaf)- The chains interact antiferromag- 
netically {Jaf)- (b) The extreme limit of the model when Jp — > 00: antiferromagnetic honeycomb 
lattice with 5=1. 

Phc. 2: Two-terminal graph used for renormalization purposes. 

Phc. 4: Neel temperature T/v vs A phase diagram for different values of the C\ = J' af /Jaf ratio: 
1.0 (1); 0.5 (2); 0.3 (3). The region above (below) the critical line represents disordered (ordered) 
phase. The dotted lines are guides-to-the-eyes. 



Phc. 5: Phase diagram Tv vs A for small T/v near A c : C\ = 1.0 and A c ~ 0.46 (a), C\ = 0.5 and 
A c ~ 0.62 (b), Ci = 0.3 and A c ~ 0.76 (c). 



Phc. 3: Critical inverse temperature (K c ) vs C\ = j' AF / Jaf found from the RG recursion relations 
for different anisotropy: A = 1.0 (1), A = 0.8 (2), A = 0.6 (3). 



